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Abstract. 

We evaluate the information capacities of a lossy bosonic channel with correlated noise. 
The model generalizes the one recently discussed in [Phys. Rev. A 77, 052324 (2008)], where 
memory effects come from the interaction with correlated environments. Environmental 
correlations are quantified by a multimode squeezing parameter, which vanishes in the 
memoryless limit. We show that a global encoding/decoding scheme, which involves input 
entangled states among different channel uses, is always preferable with respect to a local 
one in the presence of memory. Moreover, in a certain range of the parameters, we provide 
an analytical expression for the classical capacity of the channel showing that a global 
encoding/decoding scheme allows to attain it. All the results can be applied to a broad class of 
bosonic Gaussian channels. 

PACS numbers: 03.67.Hk, 03.65.Yz, 89.70.-a 

1. Introduction 

One of the main tasks of quantum information theory is the evaluation of the capacities 
of quantum channels for the transmission of classical or quantum information. Recently, 
a growing attention has been devoted to the study of quantum channels with memory. 
Coding theorems were provided for a subset of memory channels (Q, the so-called 'forgetful 
channels'. One can distinguishes the cases in which the output at the kth use of the channel 
is influenced by the input at the k'th use, with k' < k, as the models studied in [2]; and 
those in which memory effects come from correlations among subsequent channel uses, as 
the ones considered in |f3l|4l[51[6l. Here we consider the second case, which is also referred 
to as 'channel with correlated noise'. A correlation-free channel can be considered as an 
ideal limit since correlations are unavoidable in physical realizations. Another motivation 
for studying channels with correlated noise is the possibility of enhancing the information 
capacities. There are indeed evidences of the possibility of amplifying the classical capacity 
in both the cases of discrete (3]| and continuous flU variables quantum channels. 

Here we consider a model of bosonic Gaussian channel in which memory effects come 
from the interaction with a bosonic Gaussian environment. The model is a generalization of 
the one discussed in [5] and it belongs to a family of channels presented in 0. Even though 
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each channel belonging to this family is unitary equivalent to a memoryless one (in the sense 
specified in [6]), the presence of energy constraints can break the unitary symmetry, leaving 
the problem of capacities evaluation open. An instance of a channel belonging to that family is 
obtained by specifying the state of the environment. Here we consider a multimode squeezed 
thermal state, determined by two parameters. The first parameter expresses the degree of 
squeezing, which in turn determines the amount of correlations in the channel; the second 
one is a temperature parameter expressing the mixedness of the state. It is clear that at zero 
temperature the correlations in the multimode squeezed state are quantum, on the other hand 
above a certain temperature the states becomes separable and the correlations are classical. 

The choice of a Gaussian state for the environment makes in turns the channel Gaussian. 
In this way, using flTJ [8l IH [TOj, we are able to evaluate, analytically or numerically, 
the classical and quantum capacities of the memory channel. To emphasize the role of 
correlations, we compare two different scenarios for encoding and decoding classical and 
quantum information: in the first one, which we refer to as the global scenario, we allow 
preparation of states at the input field which are entangled among different channel uses; 
in the second one, called the local scenario, we only allow preparation of simply separable 
states (i.e. uncorrected) at the input field, moreover we do not allow the receiver to access the 
correlations among the output modes. 

The paper develops along the following lines. In section [2] we present the model 
and define the global and the local encoding/decoding scenarios. In section [3] we present 
analytical and numerical results for the classical, entanglement-assisted, and quantum 
capacity. Conclusions and comments are drawn in section HI 

2. A model of lossy bosonic Gaussian channel with correlated noise 

We consider an instance of the general model for a bosonic channel with correlated 
noise presented in [6]. For any integer n, its action is defined over a set of n input 
bosonic oscillators, with canonical variables {qk,Pk}k=i,...n- A collection of ancillary modes 
{Qk, Pk}k=i,...n, which play the role of the environment, is also needed. In the following we 
refer to this set of oscillators as the input and environment local modes. All the frequencies 
are assumed to be degenerate and equal to one, together with h — 1. The integer k labels the 
sequential uses of the channel. At the A;th use, the fcth input mode is linearly mixed with the 
kth mode of the environment at a beam splitter with given transmissivity 77 (see figured])- In 
the Heisenberg picture the channel transforms the input field variables as 



l'k = Vvik + y/l -riQ k , 

p'k = \Fnvk + V 1 - v p k ■ 

A constraint on the energy is required to avoid infinite capacities. We constraint the average 
number of photons at the input field; for a given N we require: 




(2) 



k=l 
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Figure 1. A schematic picture of the model of lossy bosonic channel. Each input mode (left- 
right line), representing one use of the channel, interacts with the corresponding environment 
mode (top-bottom line) through a beam-splitter. To introduce correlations effects, environment 
modes are considered in a correlated state. 



Memory effects appear in the presence of correlations among the local modes of the 
environment. For a given integer n, a channel over the n input modes is defined. In the 
Schroedinger picture its action is 

£ {n) ( Pin ) = tr cnv (U P - m ® Penv W f ) (3) 

where pi n indicates the state of the input field, U is the unitary transformation at the n beam 
splitters, p cnv indicates the state of the environment field and tr cnv the partial trace over the 
environment variables. The channel is correlation-free if the state of the environment is simply 
separable (i.e. uncorrelated) in the basis of the local mode. 

We assume the environment to be in a Gaussian state, which in turns makes the channel 
Gaussian. Here we consider an environment covariance matrix of the following block- 
diagonal form: 

(QQ T ) (S^£> 
( QpT + PQT ) (pp t ) 



v 




(4) 



where Q := {Qi,Qi 



This is a bona fide covariance 



.QnY andP := (P 1; P 2 ,...P n 
matrix as long as the matrix Vt is symmetric and T > 0. To fix the ideas we chose a n x n 
matrix Vt of the following form: 

\ 




n 



( o 
i 







1 







\ o o 














1 



1 

0/ 



(5) 



In this way, for T = 0, we recover the model analyzed in (5]|. The physical interpretation is 
that of a multimode squeezed thermal state, that can be obtained by applying the multimode 
squeezing operator (see [11]) 



U(s) = exp 



-i- > {QkPk+i + PkQk+i) 



'2^ 

k=l 



(6) 
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to a thermal state with T average excitations per mode. Since the memory of the channel 
comes from squeezing in the environment we refer to the parameter s (or \s\) as the memory 
parameter, for s = the environment is a correlation-free thermal state. The integer n can be 
interpreted as the characteristic length of correlations. To evaluate lower bounds we need the 
eigenvalues and eigenvectors of the matrix Vt which were presented in Q. The eigenvalues 
are 

Xj = 2 cos (z^lj for 3 = !> • • • n, 0) 
where the corresponding eigenvectors have components 



as: 



v i> k = \l — rr sin ( ~^~~T ) for k = l,...n. (8) 
V n + 1 \n + l j 

Let us introduce the global variables {Qj, Pj}j=i,... n , defined from the eigenvectors of ft 

Qj '■= J2k v j,kQk, ^ 

It follows that the environment covariance matrix is diagonal in this basis; it can be written as 
a direct sum 

n 

V = ®Vi (10) 

J'=l 

of single-mode covariance matrices of the form: 

)2\ . • I'iQr 




^ = ( { sMa 1 ) ( " fa ' i = ir+;n "„ i. ait 

with sj := sXj. Hence, moving to the global variables, the state of the environment is the 
direct product of n modes, each being in a squeezed thermal state, with squeezing parameter 
Sj and T thermal photons. 



2.1. Global versus local scenario 

In the following sections we estimate the quantum and classical capacities for a given value 
of the correlation length n. To emphasize the role of correlations we compare two different 
encoding/decoding scenarios. 

The first scenario is a global one, in the sense that it involves preparation of the input field 
in states which are (in general) entangled among different uses of the channel (i.e. among the 
local modes). The second scenario is a local one, involving preparation of the input field in 
states which are simply separable among the channel uses; moreover we do not allow the 
receiver to access the correlations among different local modes at the output field. 

As to the global scenario, we introduce the following global variables at the input field 

Qj '■= J2k v j,k(lk, 

Pj ■= T,k v J,kPk, 
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(notice that this set of variables is 'parallel' to the environment global variables defined in ©) 
and consider states of the input field which are factorized in the basis of these global modes. 
These states are in general entangled among channel uses (i.e. among local modes). In terms 
of the global variables, the channel factorizes as 

£^=(g)£f, (13) 

j 

where the channel 2,p acts on the jth global mode of the input field. In the Heisenberg picture 



£j transforms the jth global input variables as: 



q'j = VvQj + V^v Qj, (14) 

Furthermore, due to the form of the transformation dT2|) (a linear passive one, see 



Appendix A), the energy constraint is preserved in the basis of global variables: 



-E 



1 



J*+P 2 j )in<N+-. (15) 
3=1 

It is worth noticing that the channel £^ is a lossy bosonic channel in which one (global) mode 
is mixed with an environment mode which is squeezed (it is described by the covariance 
matrix (flTT)). Hence the n-use channel is unitary equivalent to a correlation-free channel, 
which is the product of n single-mode channels. 

Concerning the local scenario let us say that, as long as the correlations among the local 
modes at the output field are neglected, the environment can be effectively described by a 
state which is factorized in the local modes. For each k, by integrating the environment 
Wigner function over the local variables {Qh, Ph] for h ^ k,we obtain the Wigner function 
of the kth local mode of the environment. The corresponding state is thermal-like, the average 
number of photon can be computed from ©, yielding: 



T.s(k) =[T + \ 



E 

.3=1 



Vj, k \ 2 e s > 



(16) 



(2T + 1) [J2%\ Kfc| 2 cosh a A - \ if n is even, 

(2T + 1) f 2^=i \vj,k\ c oshs j) + W(n+i)/2,k\ ^ if n is odd. 



By the symmetries of Vj^ and Sj, this can be rewritten as follows: 



Hence in the local scenario the state of the environment can be substituted with a thermal state 
with ^-dependent temperature. Notice that, as one can expect, the local temperature T cS (k) 
monotonically increases with \s\. 



3. Evaluating capacities 

This section is devoted to the evaluation of classical and quantum capacities of the channel 
£H This requires the constrained optimization of several entropic quantities, as the Holevo 
information, the quantum mutual information, the coherent information. 
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Let us recall that the von Neumann entropy S(p) of an mode Gaussian state p can be 
computed as follows: 

n 

S(p) = - 1/2) (18) 



a=l 



where 



g(x) := {x + 1) log 2 (x + 1) — x log 2 x (19) 



and v a are the n symplectic invariants, i.e. the symplectic eigenvalues of the covariance matrix 
(see e.g. lfT2lO . In the case of a single mode with covariance matrix a the only symplectic 
invariant is det(cr), and the symplectic eigenvalue is v = a/ det(a). It is worth remarking that 
both in the global and local scenarios introduced above the evaluation of the entropy of the 
n-mode fields reduces to the single-mode case. 

As to the evaluation of the classical capacity of a quantum channel £, one is led to the 
Holevo information, defined as 

X(£, { Pa , d Pa }) := S(2(p)) - J d Pa S(£(p a )), (20) 

where p a denotes the quantum state encoding the classical variable a, with probability density 
dp a , and p is the ensemble state, p = J dp a p a . 

The evaluation of quantum capacity and entanglement assisted classical capacity involves 
the coherent information, defined as 

J(2,p):=S(2(p))-S(2,p). (21) 

The quantity denoted S(£, p) is the entropy exchange, defined as 

S(£,p):=S(£®3(p)), (22) 

where p is a purification of the state p involving an ancillary system, and 3 is the identical 
channel acting on the ancilla. 

It is worth remarking that Gaussian encoding it is known to be optimal for classical jHHH 
and quantum capacities [10] in the memory less case. Motivated by these results we are going 
to estimate the capacities of the memory channel using Gaussian encoding. 

Concerning the global scenario, the maximization of the entropic functions is performed 
over a set of Gaussian states which are simply separable with respect to the global modes. 
Notice that these states are in general entangled in the local basis, i.e. entangled among 
channel uses. In particular, we consider covariance matrices of the form 

n 

^ = 0^, (23) 

3=1 

where aj is the covariance matrix of the jth global input mode. For real Tj and tj > 0, the jth 
covariance matrix is chosen as follows: 
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Under the action of the channel £j- this matrix is mapped into the covariance matrix 

a] = Wo + (1-7?)^. (25) 
The energy constraint can be written in terms of the input covariance matrix as 

1 n 1 

-$>(a J )<iV+- (26) 



2n 



that is 



3=1 V 7 



coshr^iV+i. (27) 



It can be useful to write the energy constraints in two steps, namely 

t i + Ijcoshr, = Ar i + I (28) 

1 " 

-J2N,<N. (29) 

3=1 

Concerning the local scenario, as noticed above, the channel reduces to a correlation-free 
channel with thermal environment, with a A;-dependent effective temperature. For this kind 
of channel, expressions for the (one-shot) classical and quantum capacities are available in 
literature d HQl and will be used as a term of comparison. 



3.1. Holevo information 

First, let us compute an upper bound for the classical capacity of the memory channel. This 
can be obtained by the maximal output entropy. For given n, we have 

C < -sa V S{£F>(p)) < sup =: C>, (30) 

n n j= i 

where the superior is over all input states p (pj) satisfying the energy constraints, which is 
reached by Gaussian input states. The maximum output entropy is reached by input states as 
in (124)) . The contribution of the jth global mode to the output entropy reads 

S(2f\ Pj )) = g (^/det (<# - 1/2) , (31) 

where 

det (0$) = {rj{tj + l/2)e r ; + (1 - rj) (T + 1/2) e*') 

x (n(tj + l/2)e" r ' +(l-r])(T+ l/2)e- s >) . 

For sufficiently small values of s, it is possible to write an explicit solution. The maximum 
output entropy is reached in correspondence of the optimal values satisfying 

(V + 5) = *T + 5 =F ^ (T + i) sinh .„ (32) 

with 

iV° pt = iV - ^ (t + cosh Sj + l^L (M(s, T) + ~) , (33) 
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M(s, T) :=\{t + cosh " \ (34) 

The range in which these values are optimal is determined by the relations t° pt > 0, namely: 

^) 2 -(^) 2 H)^S>i (35) 
In this range the upper bound can be computed analytically, yielding 

C> = g[ V N+(l-r ] )M(s,T)]. (36) 

Let us now compute a lower bound for the classical capacity. Motivated by the fact that 
it is optimal for the memory less channel JH1 we consider encoding in displaced states: 

Qa = V(a)gV ] (a). (37) 

Here q denotes a Gaussian seed state, not necessarily the vacuum, of the n mode input field. 
As to the global scenario, its covariance matrix is assumed to be as in equations (1231) . (l24j) . 
The displacement operator can be decomposed in the basis of global modes as 

n 

Tf(a) = ®T> j {a j ), (38) 

j"=i 

where atj = (y q j + iy p> j)/\^2 is the displacement amplitude of the jth global input mode, 
and a := (a±, . . . a n ). The classical noise is assumed to be Gaussian with zero mean and 
covariance matrix Y of the following form: 

n 

Y = (BY S , (39) 

i=i 

where, for c gj , c p j > 0, 

Y ( (vh) (Vj^) ) : ( ^ ) . (40) 

At the jth global mode the classical noise induces an ensemble state described by the 
covariance matrix 

aj := <Tj + Yj, (41) 
that, under the action of the channel £^ , is mapped into 

a f j = ri^ + {\-ri)V j . (42) 

In our setting, the lower bound for the classical capacity is computed by maximizing the 
Holevo information 

n n 

X = E * = E 9 (^det(^) - 1/2) - g (^det(o^) - 1/2) (43) 
j=i 3=1 
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over the An parameters {tj, rj, c q j, c p j} under the energy constraints 

2 ( C U + c p,i) + (*i + V2) cosh = Nj + - (44) 
1 " 

-ViV,<iV. (45) 
i=i 

For given values of N, rj E (0, 1), T and for sufficiently small values of s we can write 
an explicit solution. The maximum Holevo information is reached in correspondence of the 
following optimal values of the parameters: t opt = 0, r opt = Sj, 

<% = n ^ + 1~t- 1 t-{ t+ 1) sinh s > (46) 

and 

«? = *r + \ - e ^r + l -^r ( T + 1) sinh a * (47) 

where iV° p is as in equation d33l) . 

The range of s for which these values are optimal is defined by the conditions c° pt > 0, 
c° pt > 0. In that range, we are able to provide the following analytical lower bound for the 
classical capacity per channel use: 

C< = g [ V N + (1 - r))M(s, T)} -g[(l- r/)T] . (48) 

For T = 0, this improves the lower bound computed in Q. From the lower bound we can 
deduce that, besides the trivial cases 7/(1 — 7/) =0, the classical capacity monotonically 
increases with \s\. Moreover, it is worth noticing that, at T = 0, the lower bound (|48T) 
coincides with the upper bound (1361) . Hence, in the intersection of their ranges of validity 
the analytical upper and lower bounds are strict and the expression in (l48l is the capacity 
of the memory channel at zero environment temperature. It is easy to see that the range of 
validity of the upper bound (1351) contains the one of the lower bound, thus the lower bound 
in (1481) is strict in the whole range of its validity at T = 0. Notice that the results can be 
extended to the limit n — > oo as in O. For higher values of \s\ one can look for a numerical 
solution. Figure [2] shows the analytical and numerical lower bound as function of the memory 
parameter for several values of 7/ and T. 

It is interesting to consider the limit \ s\ — > oo corresponding to 'infinite correlations'. Let 
us consider the term Xj m (1431) coming from the contribution of the jth global mode. Without 
loss of generality, we can assume Sj > 0. In the limit sj » lwe can write the following 
asymptotic expressions: 

det(^) ~ 7/(1 — n) (t+ ^ L + ^ e"- + 0(e~% (49) 



det(^) ~ 7/(1 - 7]) (t 



e s i + 0(e- Sj ). (50) 



By noticing that 

lim [g(x) - (log 2 x - log 2 e)] = 0, (51) 




Figure 2. The plots show the lower bounds for the classical capacity, for n = 10, as function 
of the memory parameter |s|. In (a) at T = 0, where the lower bound is strict, for different 
values of rj, from bottom to top r) varies from 0.1 to 0.9 by steps of 0.2. In (b) at r\ = 0.9 
for different values of T, from top to bottom T varies from to 5 by steps of 1. The solid 
lines refer to the global scenario, the tick ones refer to the analytical solution in the region in 
which it is available. The lower bounds for the local scenario are plotted in dashed lines. The 
maximum average number of excitations per mode in the input field is N = 8. 



we obtain the following asymptotic expression for the jth term in the Holevo information: 

*~> = !S. x ' = k los >( 1 + l?£j2)- <52) 
For given value of Nj, from the last expression one obtains that the maximum of the mutual 
information is reached for t opt = 0, r° pt = \n(2Nj + 1), c° pt = 0, and c op - = sinhr° pt , 
yielding the following value for the jth contribution to the classical capacity: 

C<.= max v°° = log 2 (2iV, + l). (53) 

{tj ,rj ,C q J ,c Pt j } 

Summing over j we obtain the following expression for the capacity per channel use: 

c < = [ !og 2 ( 2N + !) if n is even > (54) 

00 \ 2=1 log 2 (2N + l) + ±{g [rjN + (1 - V )T] -g[(l- V )T]} if n is odd. 

The presence of an extra term for odd n comes from the contribution of the Sj = term and it 
leads to the oscillations of the Holevo information with the number of uses already observed 
in (5). However, the relative amplitude of these oscillations becomes negligible as the number 
of uses increases. Interestingly enough, in the limit of perfect memory the maximal Holevo 
information is determined solely by the value of N, i.e. by the energy constraints. The 
asymptotic lower bound can be reached by homodyne detection (see lfl"3lO . 

To conclude this section, let us mention that a lower bound concerning the local scenario 
can be obtained by the following expression (see [6]): 



f n 1 " ] 

- max V 9 [T,iV fc + (1 - ^^(fc)] - 9 [(1 - ?,)r eff (fc)] | - Y^N k = N\. (55) 
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This lower bound saturates the channel capacity for T eS = 0, see (SJ, which is obtained for 
s = 0, T = 0. This bound is plotted in figure [2] together with the lower bound computed in 
the global scenario. 



3.2. Coherent information, quantum mutual information 

The problem of evaluating the quantum capacity is greatly simplified by the fact that the 
channel £( re ) is degradable for 77 e [1/2, 1] and anti-degradable for 77 6 [0, 1/2). It follows 
that the coherent information is additive for 77 G [1/2, 1] and the quantum capacity vanishes 
for rj E [0, l/2[ (see ||T4l IT5l ITOll X It is easy to recognize that this property is shared by all the 
Gaussian memory channels of the kind presented in [6|. For the same reason, in the global 
scenario, the n-mode channel reduces to the single-mode case as 



sup J(£W,p) = sup^J(^ 1) ,p,), 
p {Pj} ,-=1 



(56) 



with the proper energy constraint. 

For the quantum capacity, and for the entanglement-assisted classical capacity, we need 
to evaluate the entropy exchange of the channel £j. . It follows from IfTOll . and from (71 lU that 
it is sufficient to consider Gaussian input states. Numerical analysis shows that the choice of 
input state with covariance matrix of the form (1241) is optimal. The input state at the jth global 
mode, with covariance matrix as in (f24)) . can be purified into a two mode Gaussian state with 
covariance matrix: 

\ 

(57) 










X 3 










h 





-Xj 




• ? J 





h 





1. 





— Xj 








where 



"3 •- 



(t,- + l/2)e r ', bj :={tj + l/2)e 



x 



3 ■' 



a 3 b 3 



1/4. 



(58) 



The action of the channel £j ® 3 leads to the output covariance matrix: 



( rjaj + (1 - 


y/fjXj 



r]bj 



\ 







where 



"3 ■- 



U 3,± 





+ (1 



-y/rjXj 
l/2)e- 



V)dj 













(59) 



{T + l/2)e s \ dj = (T 
The symplectic eigenvalues of the covariance matrix in ( 1591) are: 

1 

V2 

where 



4det 



L 3 •' 



det (Aj) + det {Bj) + 2 det (Cj). 



(60) 



(61) 



(62) 
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(a) (b) 




S S 



Figure 3. The plots show the numerically evaluated quantum capacity, for n = 10, as function 
of the memory parameter \s\. In (a) at T = for different values of 77, from bottom to top 77 
varies from 0.6 to 0.9 by steps of 0.1. In (b) sir] = 0.9 for different values of T, from top to 
bottom T varies from to 1.5 by steps of 0.5. The solid lines refer to the global scenario, the 
dashed lines to the local one. The maximum average number of excitations per mode in the 
input field is N = 8. 



Hence, the contribution of the jth global mode to the coherent information reads 

Jj = 9 (y/detty) - 1/2) - g(u jt + - 1/2) - g{v^ - 1/2). (63) 

The constrained maximization of the total coherent information gives the quantum capacity 
of the memory channel per channel use: 




Q = - max { } } , (64) 

where the maximum is over the parameters {rj, tj} under the energy constraints (ITTT) . The 
results of numerical optimization are plotted in figure [3] 

A simple expression can be written in the limit of infinite correlations \s\ ^> 1. For 
sj > 0, we can write the following asymptotic expressions for the symplectic eigenvalues: 



7/1 ~ 
— 



'r/(l - 77) ( tj + I ) ( T + I ) e s i~ r i + 0(e-^ /2 ) (65) 



2/ V 2 



and 



Vj- ^ 2 +0{e~ s J /2 ). (66) 

Analogous expressions can be obtained for Sj < 0. Taking in account the asymptotic 
expression in (|49| ), it follows that the coherent information vanishes in the limit \sj\ — > 00: 

Jooj = lim Jj = 0. (67) 

s— too 

Hence, we can write the following expression in the limit of infinite correlations: 

if n is even, 

voo = \ 5 . f . AA (68) 
- if n is odd. 
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The finite term 



, D + N'-N-1 \ f D - N ' + N -r 
8 = g{N )-g[ - g[ , (69) 



where 



and 



N' : = r)N + (1 - rj)T (70) 



D := y/(N + N'+ l) 2 - 4r)N(N + 1), (71) 

comes from the contribution of the global mode with Sj = (see 0); however, this 
contribution becomes negligible if n >> 1 . 

The entanglement-assisted classical capacity is obtained maximizing the quantum mutual 
information: 

C e = - max If" I, I (72) 

where 

ij = g(tj) + Jy (73) 

The results of numerical maximization are plotted in figure |H In the limit of infinite memory, 
and for Sj ^ 0, the contribution of the channel £^ to the quantum mutual information is 

/oo = lim Ij = g(t 3 ). (74) 

\Sj |— >oo 

Notice that this asymptotic expression is independent of the transmissivity r] and the 
temperature parameter T. Summing over j we obtain: 



g(N) if n is even, 
#(iV) + - if n is odd. 



As for the quantum capacity the term 5 comes from the contribution of the global mode with 
Sj = 0, this contribution becomes negligible if n ^> 1. 

Concerning the local scenario, figures [310] show in dashed lines the quantum and assisted 
capacity computed applying the formulas in . 



4. Conclusion and comments 



We have presented analytical and numerical results for the capacities of a lossy bosonic 
Gaussian channel with correlated noise. To emphasize the role of correlations, we have 
compared two different scenarios. The global one allows preparation of states at the input 
field which are entangled among different channel uses. 

For our channel model we have shown that the global scenario is optimal in the presence 
of memory. In particular, we have shown that, in a certain range of the parameters, it allows 
to enhance the classical capacity over the memoryless (correlation-free) channel. The optimal 
seed state of equation (1371) turns out to be entangled as shown by figure |5] where the von 
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Figure 4. The plots show the numerically evaluated entanglement-assisted classical capacity, 
for n = 10, as function of the memory parameter \ s\. In (a) at T = for different values of 
f], from bottom to top r\ varies from 0.1 to 0.9 by steps of 0.2. In (b) at rj = 0.9 for different 
values of T, from top to bottom T varies from to 6 by steps of 1. The solid lines refer to the 
global scenario, the dashed lines to the local one. The maximum average number of excitations 
per mode in the input field is N = 8. 



4 




Figure 5. The plot shows the von Neumann entropy of the single-mode reduced state, obtained 
from the optimal seed state (see (13710. averaged over all the 1 : (n — 1) partitions. This is for 
n = 10, at T = and r\ = 0.9. The tick line refers to the analytical solution. The maximum 
average number of excitations per mode in the input field is N = 8. 



Neumann entropy of the reduced state, averaged over all the 1 : (n — 1) partitions, is plotted. 
The global scenario also allows to enhance the entanglement-assisted classical capacity, at 
least for 77 < 0.5. Moreover, it slows down the decrement of the quantum capacity, being the 
latter a decreasing function of the memory parameter. It is worth noticing that all the results 
can be generalized to a broad class of bosonic Gaussian channels with correlated noise. This 
channels are those defined by an environment covariance matrix that can be diagonalized by 
a transformation which is symplectic and orthogonal (see Appendix A). 

Finally, we comment on the role of the environment temperature. From the analytical 
and numerical results we can deduce that increasing the temperature of the environment is 
qualitatively equivalent to decreasing the beam splitter transmissivity. Hence, by increasing 
the environment temperature more noise is injected in the channel without any qualitative 
change in the behavior of its capacities. It is also worth noticing that, for fixed nonvanishing 
value of the squeezing parameter s, by increasing the temperature the environment state makes 
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(b) 




10 




Figure 6. For n = 2, the contour plots show: (a) the lower bound for the classical 
capacity, (b) the von Neumann entropy of the reduced state of the optimal seed state, (c) 
the quantum capacity and (d) the entanglement-assisted classical capacity, as function of the 
parameters determining the state of the two-mode environment: the memory parameter and 
the temperature parameter. The value of the transmissivity is r} = 0.9, the maximum average 
number of excitations per mode in the input field is TV = 8. The black line indicates the 
boundary between the region in which the state of the environment is separable (on the left) 
and entangled (on the right). 



a transition from entangled to separable. However, from the point of view of the channel 
capacities we do not find any evidence of this transition. In particular, all the (analytical 
and numerical) results are smooth functions of the environment parameters, moreover no 
qualitative difference is found in the pattern of the capacities at the transition form classical 
to quantum correlations. As an illustrative example we present in figure [6] some plots for the 
case of two channel uses. 
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Appendix A. Capacities of a broad class of bosonic Gaussian channels with correlated 
noise 

In the global scenario the channel 2> n > is unitary equivalent to a correlation-free channel 
which is the product of n single-mode channels. This equivalence was already discussed 
in j6l, however here the unitary equivalence preserves the form of the energy constraints. 
This property belongs to a large class of Gaussian memory channels. All the qualitative 
features regarding capacities are shared by all the channels belonging to this class. These 
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channels are those defined by an environment covariance matrix which is diagonalized 
by an orthogonal transformation which is also symplectic (in optics this is called passive 
transformation). Let us recall that the action of such a transformation on the phase space 
coordinates (Qi, Q2, ■ ■ ■ Qm Pi-, -^2? • • • P n ) J is by a matrix of the form (see e.g. lfT2lO 



O = „ ^ , (A.1) 




with 

XX T + YY T = I, (A.2) 
XY T - YX T = O. (A.3) 

It follows that the global scenario discussed here for the model defined by the environment 
covariance matrix in © can be equally introduced for all the covariance matrices of the 
following form 

' XD Q X T + YD P Y T YD P X T - XD Q Y T 
XD P Y T - YD Q X T XD P X T + YD Q Y T 

where the diagonal matrices Dq, Dp satisfy the Heisenberg principle, namely 

DqD p > I (A.5) 

The class of covariance matrices of the form (IA.4I) contains all the pure state covariance 
matrices and all the mixed states which are obtained by applying a squeezing transformation 



to a thermal state. For these states the global scenario can be defined as in section 12. It all 
the results concerning the global scenario can be straightforwardly extended, including its 
optimality to achieve the classical and quantum capacity. 
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